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Abstract. In the framework of the theory of open systems based on completely 
positive quantum dynamical semigroups, we give a description of the dynamics of 
entanglement for a system consisting of two uncoupled harmonic oscillators interacting 
with a thermal environment. Using Peres-Simon necessary and sufficient criterion for 
separability of two-mode Gaussian states, we describe the evolution of entanglement 
d . in terms of the covariance matrix for a Gaussian input state. For some values of the 

^ | temperature of environment, the state keeps for all times its initial type: separable or 

entangled. In other cases, entanglement generation, entanglement sudden death or a 
repeated collapse and revival of entanglement take place. We analyze also the time 
evolution of the logarithmic negativity, which characterizes the degree of entanglement 
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1. Introduction 

In recent years there is an increasing interest in using continuous variable entangled 
states in applications of quantum information processing and communication [lj. A 
full characterization of the nonclassical properties of entangled states of continuous 
variable systems exists, at present, only for the class of Gaussian states. In this special 
case there exist necessary and sufficient criteria of entanglement [2j [3] and quantitative 
entanglement measures [HE]. The quantum information processing tasks are difficult 
to implement, due to the fact that any realistic quantum system is not isolated and 
it always interacts with its environment. Quantum coherence and entanglement of 
quantum systems are inevitably influenced during their interaction with the external 
environment. As a result of the irreversible and uncontrollable phenomenon of quantum 
decoherence, the purity and entanglement of quantum states are in most cases degraded. 
Therefore in order to describe realistically quantum information processes it is necessary 
to take decoherence and dissipation into consideration. Decoherence and dynamics 
of quantum entanglement in continuous variable open systems have been intensively 
studied in the last years POmEJUniiminjES]. 
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When two systems are immersed in an environment, then, in addition to and 
at the same time with the quantum decoherence phenomenon, the environment can 
also generate a quantum entanglement of the two systems and therefore an additional 
mechanism to correlate them [91 [HI US] . Paz and Roncaglia [161 E] also analyzed the 
entanglement properties of two oscillators in a common environment by using the exact 
master equation for quantum Brownian motion and showed that the entanglement can 
undergo three qualitatively different dynamical phases: sudden death, sudden death and 
revival, and no sudden death of entanglement. 

In this paper we study, in the framework of the theory of open systems based on 
completely positive quantum dynamical semigroups, the dynamics of the continuous 
variable entanglement of two identical harmonic oscillators coupled to a common 
thermal environment. We are interested in discussing the correlation effect of the 
environment, therefore we assume that the two oscillators are uncoupled, i.e. they do 
not interact directly. The initial state of the subsystem is taken of Gaussian form and 
the evolution under the quantum dynamical semigroup assures the preservation in time 
of the Gaussian form of the state. In section 2 we write the Markovian master equation 
in the Heisenberg representation for two uncoupled harmonic oscillators interacting 
with a general environment and the evolution equation for the covariance matrix of 
the considered subsystem. By using the Peres-Simon criterion for separability of two- 
mode Gaussian states (21 HE] , we investigate in section 3 the dynamics of entanglement 
for this system. We show that for certain values of the environment temperature, the 
state keeps for all times its initial type: separable or entangled. For other values of 
the temperature, entanglement generation, entanglement sudden death or a repeated 
collapse and revival of entanglement take place. We analyze also the time evolution 
of the logarithmic negativity, which characterizes the degree of entanglement of the 
quantum state. A summary is given in section 4. 

2. Equations of motion for two harmonic oscillators 

We study the dynamics of the subsystem composed of two identical non-interacting 
oscillators in weak interaction with a thermal environment. In the axiomatic formalism 
based on completely positive quantum dynamical semigroups, the irreversible time 
evolution of an open system is described by the following general quantum Markovian 
master equation for an operator A in the Heisenberg representation (f denotes Hermitian 
conjugation) [El [20]: 

d 4 = Uh,A] + ±.J2(V?[A,VA + [Vj, A]Vj). (1) 

3 

Here, H denotes the Hamiltonian of the open system and the operators Vj, Vj , defined 
on the Hilbert space of H, represent the interaction of the open system with the 
environment. 

We are interested in the set of Gaussian states, therefore we introduce such quantum 
dynamical semigroups that preserve this set during time evolution of the system. 
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Consequently H is taken to be a polynomial of second degree in the coordinates x, y 
and momenta p x ,p y of the oscillators and Vj,Vj are taken polynomials of first degree 
in these canonical observables. Then in the linear space spanned by coordinates and 
momenta there exist only four linearly independent operators ^=1.2,3,4 [2 1 j : 

Vj = a xj p x + a yj p y + b xj x + b yj y, (2) 

where a X j, a y j, b X j, b y j are complex coefficients. The Hamiltonian H of the two uncoupled 
identical harmonic oscillators of mass m and frequency u is given by 



h = —(p 2 x + pI) + 



mu 



(x 2 + y 2 ) 



(3) 



2m yi x i y/ 2 

The fact that the evolution is given by a dynamical semigroup implies the positivity 
of the matrix formed by the scalar products of the four vectors a x ,b x ,a y , b y whose entries 

XJ ,b X j,a y j,b y j, respectively. We take this matrix of the following 
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and dissipation constant A are real 
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It follows that the principal minors of this matrix are positive or zero. From the Cauchy- 
Schwarz inequality the following relations hold for the coefficients defined in (|4]): 
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We introduce the following 4x4 bimodal covariance matrix: 
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The problem of solving the master equation for the operators in Heisenberg 
representation can be transformed into a problem of solving first-order in time, coupled 
linear differential equations for the covariance matrix elements. Namely, from ([T]) we 
obtain the following system of equations for the quantum correlations of the canonical 
observables, written in matrix form [21] (T denotes a transposed matrix): 
da it) 



dt 



Yait) + ait)Y l + 2L>, 
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The time-dependent solution of (J7J) is given by [2TJ 

a(t) = M(t)[a(0) - a(oo)]M T (t) + ct(oo), (10) 

where the matrix M(t) = exp(Yt) has to fulfill the condition lim t ^oo M(t) = 0. In order 
that this limit exists, Y must only have eigenvalues with negative real parts. The values 
at infinity are obtained from the equation 

Ya(oo) + a(oo)Y T = -2D. (11) 

3. Dynamics of two-mode continuous variable entanglement 

The characterization of the separability of continuous variable states using second-order 
moments of quadrature operators was given in Refs. [21 E]- A Gaussian state is separable 
if and only if the partial transpose of its density matrix is non-negative [necessary and 
sufficient positive partial transpose (PPT) criterion]. A two-mode Gaussian state is 
entirely specified by its covariance matrix ([6]), which is a real, symmetric and positive 
matrix with the block structure 

AC' 
C T B 



= I m d I ' (12) 



where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric 
covariance matrices for the individual one-mode states, while the matrix C contains 
the cross-correlations between modes. Simon [2] derived a PPT criterion for 
bipartite Gaussian continuous variable states: the necessary and sufficient criterion for 
separability is S(t) > 0, where 

S(t) = det A det B + Q - | det C\\ 

-Tt[AJCJBJC t J] - - (det A + det B) (13) 
and J is the 2x2 symplectic matrix 

Since the two oscillators are identical, it is natural to consider environments for which 
D xx = D yy , D xpx = D yPy , D PxPx = D PyPy , D xpy = D yPx . Then both unimodal covariance 
matrices are equal, A = B, and the entanglement matrix C is symmetric. 

3.1. Time evolution of entanglement and logarithmic negativity 



In order to describe the dynamics of entanglement, we use the PPT criterion [21 [18] 
according to which a state is entangled if and only if the operation of partial transposition 
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does not preserve its positivity. Concretely, we have to analyze the time evolution of the 
Simon function S(t) ffl~3j) . For a thermal environment characterized by the temperature 
T, we consider such diffusion coefficients, for which 

muD xx = — = — coth -77ZT, D xv = 0, 

m 2 ui 2 D xy = D PxPy . (15) 

This corresponds to the case when the asymptotic state is a Gibbs state [20] . 

For Gaussian states, the measures of entanglement of bipartite systems are based 
on some invariants constructed from the elements of the covariance matrix j6j Ell E2] ■ 
In order to quantify the degrees of entanglement of the infinite-dimensional bipartite 
system states of the two oscillators it is suitable to use the logarithmic negativity. For 
a Gaussian density operator, the logarithmic negativity is completely defined by the 
symplectic spectrum of the partial transpose of the covariance matrix. It is given by 
L = max{0, — log 2 2i>_}, where V- is the smallest of the two symplectic eigenvalues of 
the partial transpose a of the 2-mode covariance matrix a : 



Iv 2 - = A =F Y A 2 — 4 det a. (16) 

Here A is the symplectic invariant (seralian), given by A = det A + det B — 2 det C. 
In our model, the logarithmic negativity is calculated as 

2#) = ~bg a [4/(<Kt))], (17) 



where 



f(a(t)) = -(det A + det B) - det C 
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-detcx(f) • (It 



It determines the strength of entanglement for L(t) > 0, and if L(t) < 0, then the state 
is separable. 

In the following, we consider the behaviour of the logarithmic negativity for two 
types of the initial Gaussian state: 1) separable and 2) entangled. 

1) In Figures 1 and 2 we represent the dependence of the logarithmic negativity L(t) 
on time t and temperature T for a separable initial Gaussian state (unimodal squeezed 
state and respectively a mixed state). We notice that for relatively small values of the 
temperature T, the initial separable state (L(t) = 0) becomes entangled immediately or 
shortly after the initial moment of time t = (generation of entanglement, when the 
logarithmic negativity L(t) becomes strictly positive). For relatively large values of T, 
L(t) keeps its initial zero value and the state remains separable for all times. 

Depending on the environment temperature, we notice three situations in the 
case of a generated entanglement: a) entanglement may persist forever, including the 
asymptotic final state; b) there exist repeated collapse and revival of entanglement; 
c) the entanglement is created only for a short time, then it disappears and the state 
becomes again separable. The entanglement of the two modes can be generated from 
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Figure 1. Logarithmic negativity L versus time t and environment temperature T 
(via C = coth jpp) for A = 0.1, D xy — 0, D xPy = 0.049 and separable initial uni-modal 
squeezed state with initial correlations er^^O) = 3/4, cr Px p x (0) = 1/3, <J xp!C {Q) = 
Czy(0) = a VmPy (0) = a xPy (0) = 0. We take m = u = h=l. 




Figure 2. Same as in Figure 1, for a separable initial Gaussian mixed state with initial 
correlations a xx (0) = 1, cr pxPx (0) = 1/2, 0^(0) = a xy (0) = v PxPy (0) = v xPy (0) = 0. 

an initial separable state during the interaction with the environment only for certain 
values of mixed diffusion coefficient D xpy and dissipation constant A. 

2) The evolution of an entangled initial state is illustrated in Figures 3 and 4, 
where we represent the dependence of the logarithmic negativity L(t) on time t and 
temperature T for an entangled initial Gaussian state (unimodal squeezed state and 
respectively a mixed state). We observe that for relatively small values of T, the initial 
entangled state remains entangled almost for all times. For relatively large values of 
temperature T, at some finite moment of time, L(t) becomes zero and therefore the 
state becomes separable. This is the so-called phenomenon of entanglement sudden 
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Figure 3. Same as in Figure 1, for an entangled initial uni- modal squeezed state 
with initial correlations a xx (0) — 3/4, cr PxPx (0) = 1/3, & xpx (0) = 0, 0^(0) = 
1/2, a PxPy (0) = -1/2, a xPv (0) = 0. 




Figure 4. Same as in Figure 1, for an entangled initial Gaussian mixed state with 
initial correlations 0^(0) = 1, cr PxP:c (0) = 1/2, a xPx (0) = 0, a xy (0) = 1/2, cr PxPy (0) = 
-1/2, a xpy (0) =0. 

death. This phenomenon is in contrast to the loss of quantum coherence, which is 
usually gradual [121 123]. Depending on the values of the temperature, it is also possible 
to have a repeated collapse and revival of the entanglement. 

As expected, the logarithmic negativity has a behaviour similar to that one of the 
Simon function in what concerns the characteristics of the state of being separable or 
entangled [231 EQ E3 EE] . 
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3.2. Asymptotic entanglement 

From ( JTTj) and ( JTBT) we obtain the following elements of the asymptotic matrices 
A{oo) = B(oo) : 

777^0-^(00) = = -coth — — , 0-^(00) = (19) 

and of the entanglement matrix C(oo) : 

, , m 2 (X 2 + oo 2 )D xv + m\D xv , , 

**v(°°) = — 21/12 I 2 , — , 20 

rrr A(A^ + ur) 

o- xPy (oo) = o- yPx (oo) = , (21) 



m 2 u 2 (X 2 + uj 2 )D xv — mu 2 XD xv 
ff >-<°°> = A(A^) (22) 



Then the Simon expression f[T51) takes the following form in the limit of large times: 



4 V 2A;T y A 2 A 2 + w 2 



n 2 



v + < 23) 

For environments characterized by such coefficients that the expression S(oo) (T2"5j) is 
strictly negative, the asymptotic final state is entangled. Just to give an example, 
without altering the general feature of the system, we consider the particular case 
D xy = 0. Then, for a given temperature T, we obtain that 5(oo) < 0, i.e. the asymptotic 
final state is entangled, for the following range of values of the mixed diffusion coefficient 
D xPy : 

10 2D xn lu 

coth — - 1< - xpy < coth — + 1. 24 
2kT VA 2 + u 2 2kT K J 

We remind that, according to inequalities (JHJ), the coefficients have to fulfill also the 
constraint 

A uj 
— coth — — - 
2 2kT 

If the coefficients do not fulfil the double inequality ( 1241) . then S(oo) > and the 
asymptotic state of the considered system is separable. 
The asymptotic logarithmic negativity has the form 

co 2D, 



-coth— >D^. (25) 



L(ooj = -log 2 



coth ■ 



'X P y 



(26) 



2kT VA 2 + co 2 

It depends only on the mixed diffusion coefficient, dissipation constant and temperature, 
and does not depend on the initial Gaussian state. 
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4. Summary 

In the framework of the theory of open quantum systems based on completely 
positive quantum dynamical semigroups, we investigated the Markovian dynamics of 
the quantum entanglement for a subsystem composed of two noninteracting modes 
embedded in a thermal environment. By using the Peres-Simon necessary and sufficient 
criterion for separability of two-mode Gaussian states, we have described the evolution 
of entanglement in terms of the covariance matrix for Gaussian input states. For some 
values of diffusion and dissipation coefficients and of environment temperature, the state 
keeps for all times its initial type: separable or entangled. In other cases, entanglement 
generation or entanglement suppression (entanglement sudden death) take place or even 
one can notice repeated collapse and revival of entanglement. The dynamics of the 
quantum entanglement is sensitive to the initial states and the parameters characterizing 
the environment (diffusion and dissipation coefficients and temperature). We have also 
shown that, independent of the type of the initial state - separable or entangled, for 
certain values of temperature, the initial state evolves asymptotically to an equilibrium 
state which is entangled, while for other values of temperature the asymptotic state is 
separable. We described also the time evolution of the logarithmic negativity, which 
characterizes the degree of entanglement. For a given temperature, we determined the 
range of mixed diffusion coefficients for which the entanglement exists in the limit of 
long times. 
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